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ABSTRACT 

We present the results of a LIGO search for short-duration gravitational waves (GWs) associated 
with the 2006 March 29 SGR 1900+14 storm. A new search method is used, "stacking" the GW data 
around the times of individual soft-gamma bursts in the storm to enhance sensitivity for models in 
which multiple bursts are accompanied by GW emission. We assume that variation in the time differ- 
ence between burst electromagnetic emission and potential burst GW emission is small relative to the 
GW signal duration, and we time-align GW excess power time-frequency tilings containing individual 
burst triggers to their corresponding electromagnetic emissions. We use two GW emission models in 
our search: a fluence-weighted model and a flat (unweighted) model for the most electromagnetically 
energetic bursts. We find no evidence of GWs associated with either model. Model-dependent GW 
strain, isotropic GW emission energy -Egw, and 7 = ficw/^EM upper limits are estimated using 
a variety of assumed waveforms. The stacking method allows us to set the most stringent model- 
dependent limits on transient GW strain published to date. We find £qw upper limit estimates (at 
a nominal distance of lOkpc) of between 2 x 10 45 erg and 6 x 10 50 erg depending on waveform type. 
These limits are an order of magnitude lower than upper limits published previously for this storm 
and overlap with the range of electromagnetic energies emitted in SGR giant flares. 
Subject headings: gravitational waves - soft gamma repeaters 
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1. INTRODUCTION 

Soft gamma repeaters (SGRs) sporadically emit brief 
(~ 0.1s) intense bursts of soft gamma-rays. Three of 
the five known SGRs have produced rare "giant flare" 
events with initial bright, short (~ 0.2 s) pulses with 
peak electromagnetic (EM) luminosities between 10 44 
and 10 47 erg s _1 , placing them among the most EM lu- 
minous events in the Universe. According to the "magne- 
tar" model SGRs are galactic neutron stars with extreme 
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Fig. 1. — SGR 1900+14 storm light curve with 1ms bins in the (15-100) keV band. Bottom plot shows a detail. Burst start times are 
estimated by fitting the steeply rising burst edges; EM fluences are estimated by integrating light curve area under each burst. A 30-bin 
running average is shown in addition to the raw light curve. Solid lines are linear fits to rising edges; the boundaries of rising edges were 
found by examining the first derivatives in the neighborhoods of the peak locations. Crosses mark burst peaks and intersections of the 
rising edge fits extrapolated to a linear fit of the noise floor measured in a quiescent period of data in the 50 s BAT sequence before the 
start of the storm. The one-sigma timing uncertainty averaged over all measurements is 2.9 ms. X-axis times are relative to 2006-03-29 
02:53:09.9 UT at the Swift satellite. 



magnetic fields ~ 10 15 G (Duncan & Thompson 19921. 
Bursts may result from the interaction of the star s mag- 
netic field with its solid crust, leading to crustal deforma- 



tions and occasional catastrophic cr acking (|Thompson fc 
Schwartz ct al. 2005; Horowitz & Kadau 



Duncan| 


19951 


2009 


1 w 


ith su 



2009 1 with su bsequent excitation of nonrad i al neutron 



star /-modes ( Andcrsson & Kokkotas 1998 



de Freitas 



Pacheco 1998 ; Ioka 2001 ) and the emission of GWs ( Uwen 



2005||Horv ath 2005; [de Freitas PachecoTT998l |loka|2001|) . 

For re views, see |Mereghetti| ( |2008[ ); | Woods &; Thom pson 

Occasionally SGRs produce many soft-gamma bursts 
in a brief period of time; such intense emissions are re- 
ferred to as "storms." We present a search for short- 
duration GW signals (<0.3s) associated with m ultiple 



bursts in th e 2006 March 29 SGR 1900+14 storm ( Israel 
et al. 12008 1 using data collected by the Laser Int erferom- 
eter Gravit ational Wave Observatory (LIGO) (Abbott 
et al. 20071. The storm light curve, obtained from the 
Bur st Alert Telescope (B AT) aboard the Swift satel- 
lite ( Barthelmy et al.|2005 ), is shown in Fig.[l] It consists 
of more than 40 bursts in ~30 s, including common SGR 
bursts and some intermediate flares with durations > 
0.5 s. The total fluence for the storm event was estimated 
by the Konus-Wind tea m to be (1 — 2) x 10~ 4 erg cm -2 in 
the (20-200) keV range flGolenetskii et al.||2006|, imply- 



ing an isotropic EM energy A'em = (1 — 2) x 10 erg at a 



nominal distance to SGR 1900+1 4 of lOkpc (source l o- 
cation and distance is discussed in Kaplan et al. (2002 1). 
At the time of the storm both of the 4 km L1GO detec- 
tors (located at Hanford, WA and Livingston, LA) were 
taking science quality data. 

We attempt to improve sensitivity to multiple weak 
GW burst signals associated with the storm's multiple 
EM bursts by adding together GW signal power over 
multiple bursts. In doing so we assume particular GW 
emission models, which we describe in the next section. 
Fig. [2] illustrates the stacking procedure using the four 
most energetic bursts in the storm. 



2. METHODS 
The analysis is perform ed by 



the Stack-a-flare 

pipeline ( Kalmus et al.|2 009), which extends the method 
used in a recent LIGO search for transient GW associ- 



ated with individual SGR bursts ( Abbott et al.|2008 ) and 
relies on an excess power detection statistic ( jAnderson 
al. 2001). To "stack" N bursts in the storm, we first 
generate N excess power time-frequency tilings. These 
are 2-dimensional matrices in time and frequency gener- 
ated from the two detectors' data streams. Each tiling 
element gives an excess power estimate in the GW detec- 
tor data stream in a small period of time dt and a small 
range of frequency 5f. The time range of each tiling is 
chosen to be centered on the time of one of the target 
EM bursts in the storm. We then align these N tilings 
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along the time dimension so that times of the target EM 
bursts coincide, and perform a weighted addition. 

Stacking significantly improves sensitivity to GW emis- 
sion under a given model. However, improving detec- 
tion probability depends upon stacking according to GW 
emission models that correctly describe nature. The 
storm light curve motivated two stacking models: a flat- 
weighted model which includes the 1 1 most energetic EM 
bursts with unity weighting factors; and an EM-fluence- 
weighted model comprised of the 18 most energetic EM 
bursts. The N = 11 cutoff in the flat model is moti- 
vated by a clear separation in EM fluence of the 11 most 
energetic bursts in the storm. Including the 18 most en- 
ergetic bursts in the fluence-weighted model accounts for 
95% of the total EM fluence of the more than 40 bursts 
in the light curve. In the fluence-weighted model, time- 
frequency excess power tilings are weighted according to 
burst-integra ted BAT counts bef ore stacking. Further 
details are in Kalmus et al.| ( 2009 1 . 

To obtain estimates of the times of EM bursts in the 
storm, we measure the intersections of the rapid rising 
edges of each burst with the light curve noise floor mea- 
sured in a quiescent period of data in the 50 s BAT se- 
quence before the start of the storm (Fig. [TJ. We cor- 
rect these times for satellite-to-geocenter times-of-flight 
using the known SGR 1900+14 sky position and Swift 
ephemeris, which vary from (17.12-17.48) ms over the 
~30 s duration of the storm. The stack-a-flare analysis 
method is robust to relative timing error s smaller than 
GW signal durations ( |Kalmus et al 2009 ) . EM fluences 
are estimated by integrating detector counts under each 
burst in the light curve. We conservatively converted 
counts to fluences using the lower bound of the Konus- 
Wind total fluence range given above. 

We divide the GW data into an on-source time region, 
in which GWs associated with the storm could be ex- 
pected, and a background region, with statistically sim- 
ilar noise in which we do not expect a GW. This is done 
after applying c ategory 1 and categ ory 2 data quality 
cuts described in Abbott et al. ( 2009 ) . The on-source re- 
gion consists of 4 s of stacked data. Each 4 s region com- 
prising the stack is centered on the time of one of the 
EM bursts included in the GW emission stacking model. 
Background regions consist of 1000 s of data on either 
side of the storm. On-source and background segments 
are analyzed and stacked identically, and the stacked 
time-frequency tilings are passed through a clustering 
algorithm resulting in lists of "analysis events." Back- 
ground analysis events due to fluctuating detector noise 
are used to estimate the significance of on-source events; 



significant events, if any, are subject to vetoes (|Abbott 
etaT1[2009 l 



Using ±2 s regions around bursts in the storm accounts 
for uncertainties in the EM burst times and a possible 
systematic delay between GW and EM emission. Al- 
though GW emission in SGRs is expecte d to occur al- 
most simultaneously with the EM burst ( |loka||2001[ ), a 
common bias in trigger times shared by all bursts in the 
stacking set of <1 s can be handled with a ±2 s on-source 

region. 

As in Abbott et al. (2008), this search targets neu- 
tr on s tar funda mental mod e ringd o wns (RDs) predicted 



( |l'998p ; [Ioka| ( |2001[ ); |Andersson| ( |2o63p as well as short- 



And e rsson fc K o kkotas| fll998|; [de Freitas Pacheco 

OCF 



duration GW signals of unknown waveform. RDs are 
targeted because /-modes are the m ost efficient GW 
emitters ( Andersson fc Kokkotas|1998 1. We assume that 
given a neutron star, /-mode frequencies and damping 
timcscalcs would be similar from event to event, and that 
unknown signals would at least have similar central fre- 
quencies and durations from event to event. 
As in 



Abbott et al. 



( 2008 1 , we thus focus on two dis- 
tinct regions in the target signal time-frequency parame- 
ter space. The first region targets ~ 1 00-400 ms duration 
signals in the (1-3) k Hz band, which inclu des /-mode 
RD signals predicted inBenhar et al. (20041 for ten real- 
istic neutron star equations of state. VVe choose a search 
band of (1-3) kHz for RD searches, with a 250 ms time 
windo w which was f ound to give optimal search sensi- 
tivity ( Kalmus 20081. The second region targets ~(5- 
200) ms duration signals in the (100-1000) Hz band. The 
target durations are set by prompt SGR burst timescales 
(5 ms to 200 ms) and the target frequencies are set by 
the detector's sensitive region. We search in two bands: 
(100-200) Hz (probing the region in which the detectors 
are most sensitive) and (100-1000) Hz (for full spectral 
coverage below the RD search band) using a 125 ms time 
window. In all, we search in three frequency bands and 
two GW emission models (flat and fluence-weighted). 
This amounts to a total of six 4 s-long stacked on-source 
regions. 



We estimate loudest-event upper limits ( Brady et al. 
20041 on GW root-sum-squared strain h TS l incident at 
the detector. We can construct si mulations of impingin g 
GW with a given h rss . Following Abbott et al.| ( |2005b[ ) 

SSX! ( 1 ) 



h 2 

"TSS + 



where e.g. 



' dt 



(2) 



and h+ jX (t) are the two GW polarizations. The rela- 
tionship between the GW polarizations and the detector 
response h(t) to an impinging GW from a polar angle 
and azimuth (#, 4>) and with polarization angle i\) is: 

h(t) = F+(0, 0, $)h+(t) + F x (9, <j>, $)h x (i) (3) 

where F + (9, <j>, tji) and F x (9, 4>, i/j) are the antenna func- 
tions for the source at (9, (f>) ( |Thorne||1987 l. At the time 
of the storm, the polarization-independent RMS antenna 
response (F? + F 2 ) 1 / 2 , which indicates the average sensi- 
tivity to a given sky location, was 0.39 for LIGO Hanford 
observatory and 0.46 for the LIGO Livingston observa- 
tory. 

We can also set upper limits on the emitted isotropic 
GW emission energy .Egw at a s ource distance R asso- 
ciated with h+(t) and h x (t) via (Shapiro & Teukolsky 



1983) 



E GW = AnR 



16ttG 



((h + r + (h x y)dt. (4) 



The procedure for estimating loudest-event up per lim- 
i ts in t he individual burst s earch is detailed in lKalmusl 
( |2008| ); |Abbott et al] ( |2008[ ). In brief, the upper limit 
is computed in a trequentist framework by injecting ar- 
tificial signals into the background data and re covering 
them with the search pipeline (see for example Abbott 
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time from stack center [s] 



Fig. 2. — Individual EM bursts inform the stacking of GW data. This figure suggests the stacking procedure and explicitly shows search 
timescales. The top four plots are EM light curve time series around individual bursts beginning at 2.0 s, 16.6 s, 18.3 s, and 22.3 s in Fig. 
^ Simulated GW ringdowns in the fluence-weighted model are superposed. The bottom plot shows the EM time series simultaneously, and 
the sum of the hypothetical GW signals. The on-source region of ±2s is shaded. In the search, GW data corresponding to the EM time 
series are transformed to time-frequency power tilings before being added together and therefore there is no dependence on phase-coherence 
of GW signals in the analysis; this transformation is not illustrated. 



et al. (2005a 2008 1). An analysis event is associated with 
each injection, and compared to the loudest on-source 
analysis event. The GW strain or isotropic energy at 
90% detection efficiency is the strain or isotropic energy 
at which 90% of injections have associated events louder 
than the loudest on-source event. 

We use th e twelve waveform types described in |Abbott1 
et al. ( 2008[ ) to establish detector sensitivity and thereby 
set upper limits: linearly and circularly polarized RDs 
with r = 200 ms and frequencies in the range (1-3) kHz; 
and band- and time-limited white noise bursts (WNBs) 
with durations of 11 ms and 100 ms and frequency bands 
matched to the two lower frequency search bands. 

These waveforms are used to construct compound in- 
jections determined by the emission model. In the flat 
model, 11 GW bursts comprise a compound injection, 
each is identical, and our stated h TSS and -Eqw are for 
one such GW burst in the compound injection. In the 
fluence-weighted model, 18 GW bursts comprise a com- 
pound injection, they are weighted (in amplitude) with 
the square root of integrated counts, and our stated h ras 
and £gw are for the loudest GW burst in the com- 
pound injection. A single polarization angle is chosen 
randomly for every compound injection. In assuming 
that the bursts emitted are identical up to an amplitude 
scale factor, we implicitly assume the star's GW emission 
mechanism and symmetry axis are constant over bursts 
in the storm. 

3. RESULTS 

We find no statistically significant GW signal associ- 
ated with the SGR 1900+14 storm. The significance of 
on-source analysis events is inferred by noting the rate 
at which background analysis events of equal or greater 
loudness occur. We examined 4s stacked on-source re- 
gions in the flat and fluence-weighted models in the three 
search bands. The most significant on-source analysis 
event from these six searches was from the flat model in 
the (100-1000) Hz band and had a corresponding back- 



ground rate of 5.0 x 10 _2 Hz (1 per 20 s) in that search. 

Tableland Fig. [3jgive model-dependent loudest-event 
upper limits at 90% detection efficiency computed for 
the GW signal associated with the single loudest EM 
burst. We give strain upper limits (h®^°) and isotropic 
emission energy upper limits at a nominal SGR 1900+14 
distance of lOkpc (E$?$). We also give upper limits 
7ul = Eqw/E-em, a source-distance-independent mea- 
sure of the extent to which an energy upper limit probes 
the GW emission efficiency, calculated using a conserva- 
tive estimate of 1.0 x 10~ 4 erg cm~ 2 for the total flucncc 
of the storm to estimate fluences for individual peaks. In 
the fluence-weighted model, 7 is the same for each indi- 
vidual burst. In the flat model we report the mean value 
of 7 for the 11 bursts. 

Superscripts in Table[T|give a systematic error and un- 
certainties at 90% confidence. (Similar estimates were 
made for the E^ff but are not shown in the table.) The 
first and second superscripts account for the systematic 
error and statistical uncertainty, respectively, in the de- 
tector calibrations. The third is the statistical uncer- 
tainty from using a finite number of trials (200) in the 
Monte Carlos, esti mated with t he bootstrap method us- 
ing 200 ensembles (Ef ron| |1979[ ) . The systematic error 
and the quadrature sum of the statistical uncertainties 
are added to the final sensitivity estimates. One-sigma 
burst timing uncertainties from fits of burst rising edges 
are accounted for in the Monte Carlo simulations. Esti- 



mating uncertainties is further described in Kalrnus et al. 
p009|. 



4. DISCUSSION 

The stacked search described here extends the recent 
LIGO search for GW associated with the 2004 SGR 
1806-20 giant flare and 190 lesser events fro m SGR 1806- 
20 and SGR 1900+14 ( [Abbott et al.|[2008] ). That search 
was the first search sensitive to neutron star /-modes, 
and it set individual burst upper limits E^S^ ranging 
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TABLE 1 

Stack-a-flare SGR 1900+14 storm upper limits. 







N= 


11 Flat 








N=18 Flucnce-wei 


ghted 




Simulation type 




/i? s ° a % [10- 22 Hz~ s] 




E%°% fcrgl 

CtW L oj 


OTJL 




h? s ° s % [10- 22 Hz" i] 




E?SS fergl 




WNB 11ms 100-200 Hz 


1.3 


+0.0 +0.17 +0.0 


= 1.5 


1.9 X 10 45 


3 x 10 4 


2.1 


+0.0 +0.27 +0.094 


= 2.4 


5.0 x 10 45 


3 x 10 4 


WNB 100ms 100-200 Hz 


1.5 


+0.0 +0.19 +0.0 


= 1.7 


2.4 x 10 45 


4 x 10 4 


2.3 


+0.0 +0.30 +0.098 


= 2.6 


6.0 x 10 45 


3 x 10 4 


WNB 11ms 100-1000 Hz 


3.5 


+0.0 +0.45 +0.0 


= 3.9 


1.8 x 10 47 


3 x 10 6 


5.2 


+0.0 +0.67 +0.29 


= 5.9 


4.1 x 10 47 


2 x 10 6 


WNB 100ms 100-1000 Hz 


3.8 


+0.0 +0.50 +0.0 


= 4.3 


2.0 x 10 47 


3 x 10 6 


5.6 


+0.0 +0.73 +0.29 


= 6.3 


4.5 x 10 47 


2 x 10 6 


RDC 200ms 1090 Hz 


4.5 


+0.045 +0.59 +U.0 


= 5.2 


1.2 x 10 48 


2 x 10'' 


7.2 


+0.072 +U.93 +U.33 


= 8.2 


3.0 x 10 48 


2 X 10' 


RDC 200ms 1590 Hz 


6.4 


+0.19 +0.84 +0.0 


= 7.4 


5.1 x 10 48 


8 x 10 7 


11 


+0.33 +1.4 +0.44 


= 13 


1.5 x 10 49 


8 x 10 7 


RDC 200ms 2090 Hz 


9.3 


+0.28 +1.8 +0.41 


= 11 


2.1 x 10 49 


3 x 10 s 


14 


+0.43 +2.8 +0.72 


= 18 


4.9 x 10 49 


3 x 10 s 


RDC 200ms 2590 Hz 


11 


+0.34 +2.2 +0.32 


= 14 


4.6 x 10 49 


8 x 10 s 


17 


+0.50 +3.3 +1.0 


= 21 


1.0 x 10 50 


5 x 10 s 


RDL 200ms 1090 Hz 


9.3 


+0.U +1.2 +u.y5 


= 11 


5.3 x 10 48 


9 x 10' 


16 


+U.U +2.1 +1.H 


= 18 


1.5 x 10 49 


8 x 10' 


RDL 200ms 1590 Hz 


14 


+0.42 +1.8 +1.1 


= 17 


2.6 x 10 49 


4 x 10 s 


19 


+0.58 +2.5 +1.9 


= 23 


5.1 x 10 49 


3 x 10 8 


RDL 200ms 2090 Hz 


20 


+ 1.2 +3.9 +1.4 


= 25 


1.0 x 10 50 


2 x 10 9 


27 


+ 1.6 +5.3 +2.8 


= 34 


1.9 x 10 50 


1 x 10 9 


RDL 200ms 2590 Hz 


25 


+ 1.8 +5.0 +3.0 


= 33 


2.6 x 10 50 


4 x 10 9 


39 


+2.7 +7.7 +2.5 


= 50 


6.2 x 10 50 


3 x 10 9 







from 3 x 10 45 erg to 9 x 10 52 erg (depending on waveform 
type and detector antenna factors and noise character- 
istics at the time of the burst), b ut did not det e ct an y 
GWs. The best values of 7 UL in |Abbott et ah] Q2008| ), 
for the giant flare, were in the range 5 x 10-6 x l(j 6 
depending on waveform type. 

The upper limits obtained here are a factor of 12 more 
sensitive in energy than the SGR 1900+14 storm upper 
limits in 
in a sing 



Abbott et al. (20081, which analyzed the storm 
e ±20 s on-source region. Those previous lim- 
its already overlapped the range of EM energies seen in 
the loudest flares as well as the range of GW en ergies 
predicted by the most extreme models (Ioka 2001). The 
flat model gives isotropic energy upper limits on aver- 
age a factor of 4 lower than a reference N = 1 (non- 
stacked) scenario (with a ±2s on-source region) and a 
factor of 2 lower than the fluence-weighted model. How- 
ever, our storm 7 upper limits are still a few hundred 
times the SGR 1806-20 giant flare 7 upper limits, due to 
the tremendous EM energy released by the giant flare. 
There is very little discussion of 7 in the theory litera- 
ture with which to compare. 

The Advanced LIGO detectors promise an improve- 
ment in energy sensitivity of more than a factor of 100. 
Furthermo re, on 2008 August 22, SGR 0501+4516 was 
discovered (|Holland et al. |2008| [Barthelmy et al.||2008| 
Palmer & Uarthclmy 2008 ) and may be located only 
1.5 kpc away (H"aensler fe Chatterjee|2008| |Leahy fc As 



chenbach||1995 1. SGR 0501+4516 searches will thus gain 
an additional 2 orders of magnitude in energy and 7 
upper limits compared to SGRs at 10 kpc. A stacking 
analysis of SGR 0501+4516 bursts with Advanced LIGO 
(a gain of 4 orders of magnitude in energy sensitivity) 
could therefore reach 7 values below unity, even without 
another giant flare. 

In the future we plan to carry out stacking searches on 
isolated SGR bursts, and eventually to perform searches 
using Advanced LIGO data. Our stacked upper limits 
depend on theoretical guidance as to what weightings 
and time delays are possible, and the significance of our 
results depends on predictions of the range of -Egw an d 
7; yet all of these things are scarce. We hope that our 
continued efforts to search for GW associated with SGR 
and Anomalous X-ray Pulsar bursts encourage further 
modeling of GW emission from these intriguing objects. 
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Fig. 3. — Stack-a-flare SGR 1900+14 storm isotropic energy 
upper limit estimates at 10 kpc, for fiat and fluence-weighted emis- 
sion models. We set upper limits at characteristic points in the 
signal parameter space in order to quantify the meaning of our 
non-detection result. Uncertainties have been folded in. Verti- 
cal lines indicate boundaries of the three distinct search frequency 
bands. Crosses and circles indicate linearly and circularly polar- 
ized RDs, respectively. Triangles and squares represent 11ms and 
100 ms band- and time-limited WNBs, respectively. Symbols are 
placed at the waveform central frequency. These results reflect the 
noise curves of the detectors. 
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